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ALEPH/OPAL data on the spectral functions from hadronic r decays are used in connection with a set of Laplace 
transform sum rules (LSR) for fixing the size of the QCD vacuum condensates up to dimension 18. Our results favour the 
ones from large- N c QCD within the Minimal Hadronic Approximation (MHA) and show a violation of about a factor 2-5 
of the vacuum saturation estimate of the dimension-six to -ten condensates. We scrutinize the different determinations 
of the QCD vacuum condensates using r-decays data. After revisiting some of the existing results, we present coherent 
values of the condensates from different methods. 



1. INTRODUCTION 

Hadronic tau decays have been demonstrated [1] 
to be an efficient laboratory for testing perturba- 
tive and non-perturbative QCD. That is due both 
to the exceptional value of the tau mass situated 
at a frontier regime between perturbative and non- 
perturbative QCD and to the excellent quality of the 
ALEPH/OPAL [2,3] data. On the other, it is also 
known before the advent of QCD, that the Weinberg 
[4] and DM0 [5] sum rules are important tools for 
controlling the chiral and flavour symmetry realiza- 
tions of QCD, which are broken by light quark mass 
terms to higher order [6] and by higher dimensions 
QCD condensates [7] within the SVZ expansion [8] 1 . 
For completing our program in the vector and V + A 
channel [10-14], we probe, in this paper, the struc- 
ture of the QCD vacuum using the ALEPH/OPAL 
data on the V-A spectral functions in connection 
with a set of Laplace sum rules (LSR). We have al- 
ready initiated the analysis of the V — A channel in 
previous papers [15,16]. However, our main motiva- 
tion here is due to the recent interests on the V — A 
hadronic correlator, which can serve as an order pa- 
rameter of spontaneous chiral symmetry breaking in 
the chiral limit m q = 0. This correlator also governs 
the dynamics of the weak matrix elements of the 
electroweak Penguin-like operators [17,18]. These 
important properties require a detailed structure of 
the related QCD vacuum which can be parametrized 
by the sum of power corrections [8]. A large num- 
ber of papers on the estimates of these power cor- 
rections using different methods exist in the litera- 
ture, but with conflicting results in [2,3,19-23] and 
in [24-26,28,29]. In the following we propose a set of 
Laplace transform sum rules (LSR) which can help 
to clarify such discrepancies. We shall be concerned 
here with the V-A two-point correlator: 
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built from the left- and right-handed components of 
the local weak current: 
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Following SVZ [8], the correlator can be approxi 
mated by: 
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where 2 d = C2d{02d) is the short hand-notation of 
the the QCD non-perturbative condensates (02d) of 
dimension D = 2d and its associated perturbative 
Wilson coefficient C 2 d] q 2 = —{Q 2 > 0) is the mo- 
mentum transfer. In the chiral limit m u ,d = 0, there 
is no D — 2 term as unflavored contribution of the 
renormalon-type [30,10,11] vanishes. The spectral 
function (v — a): 
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has been measured by ALEPH and OPAL [2,3] us- 
ing r-decay data. Within a such normalization, the 
original "sacrosante" first and second Weinberg sum 
rules [4,5] read, in the chiral limit m Ul d = 0: 
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where f n = (92.4 ± 0.26) MeV is the experimental 
pion decay constant. 

2. THE LAPLACE SUM RULES (LSR) 

In order to exploit the ALEPH/OPAL [2,3] data on 
the spectral function v — a from hadronic tau-decays, 
we shall work with the LSR version of the 1st Wein- 
berg sum rule, in the chiral limit m u ,d — 0: 



1 For a review, see e.g. [9]. 
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from which on can obtain, by taking successive 
derivatives in r, the set of LSR: 
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For our purpose, we shall truncate (in order to have 
a much better comparison with the existing results) 
the series at 2d = 18-dimension condensates 2 , as- 
suming that this approximation provides a good de- 
scription of the exact expression of the two-point cor- 
relator H LR . Then, from our previous general for- 
mula, one can write the set of sum rules: 

£ 8 ~ +018 
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• Therefore, we can extract iteratively the vacuum 
condensates beginning from £$. The value of Ci 8 
obtained in this way will be inserted into £7 for de- 
termining 0i6 and so on. 

• We parametrize the spectral function by using the 
ALEPH/OPAL [2,3] data on the spectral function 
v — a from hadronic tau-decays below t c . Above t c , 
we use a QCD continuum coming from the disconti- 
nuity of the QCD diagrams. In the particular case 
of lm!l LR studied here, a such contribution vanishes 
identically, which is equivalent to cut the integral 
in Eq. (8) at t c . The appropriate values of t c has 
been studied in [20,15,24,31] by requiring that the 
1st and the 2nd Weinberg sum rules vanishes in the 
chiral limit to leading order of the OPE. Two solu- 
tions have been found: 

t c ~ (1.4 ~ 1.5) GeV 2 and (2.5 - 2.6) GeV 2 . (9) 

In [15] the lowest value of t c has been favoured due to 
the inaccuracy of the ALEPH/OPAL data which af- 
fects the highest value, though intuitively, one tends 
to favour this highest value of t c where pQCD is ex- 
pected to work better. Exluding the high-i c value 
solution, and requiring simultaneous zeros of the 1st 
and 2nd Weinberg sum rules, Ref. [15] deduces the 
accurate number: 

t c = (1.475 ± 0.015) GeV 2 . (10) 

2 The result will not depend crucially on the choice of the trun- 
cation of the series. We shall see that at the region where the 
condensates are estimated the OPE presents a good conver- 
gence. 



This choice, as emphasized in [20] coincides with the 
t c -value obtained for the MHA in the large iV c -limit, 
which follows from the duality relation [31]: 
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with g A ~ 0.5 ± 0.06. Instead in [24 26], the higher 
value of t c around 2.5 GeV 2 has been favoured. 

• In order to avoid results which strongly depend on 
these choices of t c , we only consider the above val- 
ues of t c as a guideline of our analysis. Indeed, it 
is unlikely to take t c < 1.4 GeV 2 as we will loose 
part of the p meson tails, and then most of the low- 
est ground state dynamics. Taking t c > 2.6 GeV 2 , 
the kinematic region is small and the data become 
very inaccurate. Then, they cannot provide useful 
information to the spectral function. Indeed in this 
region, the spectral function does not have a definite 
sign, for a given data point, due to the large error 
bars. 

• For an illustration, we show the analysis of 0is,i6 
and 06,4 m Fig- 1 f°r different choices of the t c - 
cut until which we use the ALEPH/OPAL data, and 
beyond which the pQCD diagram is expected to de- 
scribe the two-point correlator. The analysis of the 
other condensates present similar features. 

• The optimal results given in Table 1 correspond to 
the one at the minimum or inflexion point of r for 
different t c -values inside the range in Eq. (9). The r- 
stability criterion has been often used in the Laplace 
sum rules analysis as it signals the compromise re- 
gion where the OPE is reliable (smaller r- values) and 
where the information from the data still remains op- 
timal (larger r- values) . It is also unlikely if the result 
is strongly dependent on the choice of £ c -values as 
this signals a strong model dependence of the result 
on the form of the QCD continuum. Then, in the 
following, we shall use in connection these two sta- 
bilities criteria for extracting the optimal results 10 . 

• The error takes into account the one of the data 
and the systematics of the method due to the range 
of t c -values given in Eq. (9) and to the propagation 
of errors induced by the ones of the input conden- 
sates. We do not include some eventual statistical 
errors. 

• It is important to notice from our analysis that 
in the range of t c given in Eq. (9), the extracted 
values of the condensates do not flip sign contrary 
to the case of FESR's results given in the existing 
literature. One can attribute this feature to the role 
of the exponential weight in LSR which enhances the 
contribution of the low-energy region to the sum rule. 

• One can also notice that, for high-dimension con- 
densates, the optimal values are obtained at large 
r-values like also in the least square fit analysis of 
[22,23]. However, we have checked that, during the 
analysis of each sum rule, the high-dimension con- 

10 For more complete discussions, see e.g. see e.g. [9]. 




Fi gure 1. t in GeV "^-behaviour of the condensates 02d i n units of GeV^ for different values of t c in GeV^: 1.4 (dot-dashed), 
1.5 (dashed bold), 2.5 (continuous bold), 2.6 (continuous): a) Oi8, b) Oie, c) Oe and d) O4. 



densates remain corrections to the low-dimension 
contributions and do not break the OPE. One can 
also notice that the position of the minimum shifts 
to lower values of r for decreasing dimension conden- 
sates, as one can see in Fig. 1 for the D = 18 to the 
D — 4 condensates. These features are re-assuring 
for the reliability of the result. 

• In order to test the accuracy of our estimate, we 
have extracted from Co the known tiny value of the 
04 quark condensate contribution using as input all 
higher dimension condensates. Including radiative 
corrections, this contribution reads [8,32]: 
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where (m u + md){uu + dd) = — 2/^m 2 .. The size of 
the radiative corrections is about 35% at the r-scale 
where the optimal results are extracted. This gives 
in units of 10~ 3 GeV 4 : 



-0.44 , (13) 
in excellent agreement with our fit — 0.5 ± 0.1 given 
in Table 1 from Fig. Id). This test increases the con- 
fidence on our other predictions in Table 1 obtained 
in the same way. 



3. ALTERNATIVE 

AND e> 8 
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(D being the dimension of the condensates). Here, 
we present alternative estimates based on some com- 
binations of LSR in the chiral limit m u ( i = 0. 
• The first sum rule is chosen in such a way that 
08 disappears to leading order while higher dimen- 
sions D = 10, 12 have smaller coefficients than in the 
individual sum rules. 
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• In the second sum rule, 06 disappears and then 8 
will dominate the LSR: 
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• Therefore, we use the sum rule in Eq. (14) (resp. 
Eq. (15)) for extracting 6 (resp. 8 ). We use the 
known tiny value of D = 4 quark condensate contri- 
bution given in Eq. (13). The analysis is shown in 
Fig. 2 and the results are given in Table 1. 

• Finally, we analyze the r-like decay sum rule, which 
has the advantage to be kinematically suppressed 
near the real axis: 



Using the previous method, we have obtained from 
Eq. (8) the results in Table 1 in units of 10~ 3 GeV 15 
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from which we deduce Oq using as input O4 and the 

higher dimension condensates. 

• Our different results are summarize in Table 1. 
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Figure 2. Same as in Fig. 1 but for the improved analysis 
in section 3: a) Os from Eq. (15), b) Oe from Eq. (14) and 
c) O g from Eq. (16). 



4. COMPARISON WITH EXISTING ESTI- 
MATES 

a) Large-A^ and Minimal Hadronic Approxi- 
mation (MHA) 

• Our results agree in signs and in magnitude until 
the D = 14-dimension condensates with the ones 
in [19] obtained using large N c and the minimal 



hadronic approximation (MHA) and with its im- 
proved version including the next radial vector me- 
son p' . 

• Our result for the D — 16 condensate still agrees 
in sign with the one in [19] but our absolute value is 
lower than the one in [19] by about 2c 

b) ALEPH and OPAL estimates from r-decay 

Our results for the low D = 6, 8 condensates agree 
also quite well with the ALEPH/OPAL estimate of 
the separate V and A channels [2,3], [21], from which 
one can deduce the V-A difference. 

c) Exponential-like sum rules 

• The value of the D = 6 condensate also agrees 
within the errors with the results in [22] , but the val- 
ues of O8,io obtained in the present paper are about 
two times higher. However, our analysis differs from 
[22] who use a least square fitting procedure with 
some other forms of LSR with a different kernel. Due 
to the alternate sign of the condensate contributions 
to the OPE, the fitting procedure can be inaccurate 
as we shall see explicitly in a forthcoming example. 

d) Finite Energy Sum Rules (FESR) 

In [24-26], FESR: 
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(n=0,l,2,3), and its slight variants have been used for 
determining 06,8- However, unlike the case of LSR 
analyzed in previous sections, the results depend cru- 
cially on the choice of t c at which one extracts the 
optimal results. The two sets of i c -values correspond- 
ing to the zeros of the 1st or/and 2nd Weinberg sum 
rules are given in Eq. (9). The results from LSR arc 
consistent with the ones corresponding to value of 
t c « 1.5 GeV 2 , while instead in [24-26], the higher 
value of t c f» 2.5 GeV 2 has been favoured. As a con- 
sequence, the value of Os and other higher dimension 
condensates obtained in these works are opposite in 
signs 11 with the ones from LSR and from MHA in 
large N c . Taking the value of t c in Eq. (10), we give 
the version of the FESR results of [24,26] in Table 
1, where the slight difference is due to the different 
parametrizations of the r-decay data (neural network 
in [26]) and to the different weights introduced for 
improving the original FESR. 

e) Weighted Finite Energy Sum Rules 

This FESR-like sum rule called "pinched-weight 
FESR" (hereafter denoted wFESR) by the authors 



11 We have corrected the sign of Os in the curve a), Fig. 5 
of [25]. Therefore curve a) and b) cross at t c fa 1.3 GeV 2 
giving the value of Os in Table 1. We have also rescaled the 
normalization by a factor 2 for consistency in our comparison. 




Figure 3. i c -bchaviour in GeV 2 of different observables used in [28]: a) (t 2 c /7)j wi , b) {t\/2)J W2 , c) 6 , d) O a in units of GeV 2d , 
2d being the dimensions of the condensates. The two curves delimit the region induced by the errors of the data. They coincide in 
almost all regions except the ones above 2.4 GeV 2 where the data are inaccurate. 



[28] is an involved variant of the FESR in Eq. (17): 
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where the weight factor u) n is: 
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for n = 2, 3, 4, 5, 6, and corresponds to the so called 
maximally safe analysis. The QCD expressions of 
these sum rules are given in Eq. (24) of Ref. [28] 
which we have checked the LO terms. 
• In order to test the results, we study the t c - 
dependence of J Ul and J UJ2 as shown in Figs. 3a) 
and 3b). We include into the analysis, the known 
effect of 04, which we have also recovered in the 
previous section. We obtain a testability point (a 
compromise region between the convergence of the 
OPE (small t c ) and minimal dependence on the form 
of the QCD continuum (large t c j) around 2.0 GeV 2 , 
at which we can extract the optimal value of the con- 
densates. However, one can notice from Figs. 3 that 
contrary to FESR, the estimates are not very sen- 
sitive (change in the last digit) to the values of t c 
corresponding to the range in Eq. (9). Neglecting 
the small radiative corrections for illustration, one 
obtains in units of 10 -3 GeV 6 : 
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We insert into this expression the values in units of 
10~ 3 of 6 , 8 fitted by [28] (mean value from ALEPH 
and OPAL fit): 



6 ~ -4.9 , 
which gives: 
i(i) 
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(21) 



0^(2.15) ~ -5.7 and 0^(2.15) ~ -5.8 . (22) 

This test shows the consistency between the results 
obtained using t c stability in Eq. (20) and the least 
square fit in Eq. (21). 

• Alternatively, we can also solve the two equations 
J UJl and J U2 for extracting the two solutions 06 and 
08- We study the £ c -dependence of the results in 
Figs. 3c) and 3d). Here the stability is obtained at 
t c ~ 1.7 GeV 2 which differs from the one obtained 
previously. We may interpret this difference as due 
to the fact that we do not consider here the same 
observables as in Figs. 3a) and 3b). However, in 
order to have conservative results, we shall consider 
t c in the range (1.7 ~ 2.0) GeV 2 where the two sta- 
bilities are obtained. One can inspect that, in this 
range, the estimate is only slightly affected by the 
t c -values. The results are given in Table 1 with a 
good accuracy. We check again the consistency of 
the results by inserting these values into Eq. (20), 
which gives: 



0$ (2.15) ~ -6.7 and 0, 



(2) 



-6.1 



(23) 
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• Our test does not support the results given in [28] 
obtained from numerical fits. This may due to the 



fact that the terms in the scries have alternate signs, 
and/or where the 2nd term is a small correction of 
the 1st one, and may be difficult to extract from 
the fitting procedure . Instead, we expect that the 
new results from this method which we give in Ta- 
ble 1 obtained using stability criteria, from solving 
the two equations J LOl and J^ 2 for extracting the two 
unknown Oq and 0§ are more reliable. 

• Notice that in a large range of t c , the two estimates 
of 06 and 0% do not flip sign, which, like in the case 
of LSR, can be due to the weight factor in the spec- 
tral integral. This is not the case of the basic FESR. 

• In principle, once one knows 06 and 08, one can 
extract the other high-dimension condensates from 
the set of equations given in Eq. (24) of [28] . 0io to 
016 can be, e.g., extracted from J U3 to J Wl0 . How- 
ever, more we go to higher moments, less the accu- 
racy on the estimate is reached as the high-dimension 
terms which one wishes to extract are tiny correc- 
tions to the leading order terms, while the method 
is not accurate enough to pick up these tiny correc- 
tions. For instance at t c « 2 GeV 2 , the QCD parts 
of the sum rules normalized to the leading 6 con- 
tributions read: 

4 



dcr in N c at Q 2 = v 2 [24]: 
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The corrections are a factor 2 smaller for J U9 and 
J Wl0 . This fact may explain why relatively large cen- 
tral values of the high-dimension condensates emerge 
from this method. A tentative extraction of 0io to 
016 from J LU3 to J Wf5 shows that the independence 
present a flat stability around 1.25 GeV 2 and another 
extremum around 1.8 GeV 2 . The values obtained at 
the second point are very sensitive to the input value 
of 06 and flip sign compared to the one at the flat 
plateau for D > 14, a feature similar to 8 from 
FESR analysis [24,26]. This may indicate that the 
weight factor is less efficient for high-dimension con- 
densates. We have excluded the high-i c solution sim- 
ilarly to the FESR case, and we deduce the values in 
Table 1. 

f) Test of the factorization assumption 

• The D — 6 condensate contributions to U LR have 
been first derived in [7] using the leading order result 
of [8] for the vector and axial- vector correlators. The 
radiative corrections have been obtained in [33,34]. 
Using an anti-commuting 75 matrix and the choice 
of operator basis in [34], it reads by assuming a fac- 
torisation of the four-quark condensates : 



64 

06 = - —Tra s (uu) 



7T 



89 
48 



(25) 



Using the NDLR or/and the HV rcgularization 
scheme, the same contribution reads, to leading or- 
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where Alr ~ (4.4 ± 0.5) x 10~ 3 is of order a 2 . 
• The D = 8 four-quark condensate contributions 
have been obtained in [35] where a 1/N 2 ambiguity 
has been noticed. The D = 10 condensates have 
been obtained in [22]. Assuming factorization, one 
can write: 
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• Mq is the scale gouverning the mixed condensate 
and is equal to (0.8 ± 0.2) GeV 2 from the baryon 
sum rules [36,37], B-B* mass-splitting [38] and string 
model [39] . We shall use the value of the gluon con- 
densate (a s G 2 ) = (0.07 ± 0.01) GeV 4 from e+e" 
data [12,40]. Within the factorization assumption, 
we shall include the log-dependence of the quark con- 
densate and of a s , which give: 
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which is 1.6 x 10~ 4 GeV 6 if one uses the invariant 
quark condensate (uu) ~ — (248 MeV) 3 [41] and eval- 
uate a s (r) at t = 1.5 GcV~ 2 at which 8 has been 
extracted from the sum rule. Using these numerical 
inputs, we deduce in units of 10 _3 GeV £) (D being 
the dimension of the condensate): 
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• Comparing these values with the ones in Table 1 
and section 3, we conclude that the factorization as- 
sumption agrees in sign with these results but un- 
derestimate the absolute value of the condensates by 
a factor 2-5. This feature is similar to the case of 
the vector [14,42,12], axial- vector [43,11], baryon [37] 
sum rules and from the analysis of the V and V + A 
r-decay data [2,3,11]. From the theoretical point of 
view, the factorization assumption is only consistent 
with the renormalization of operators to leading or- 
der in 1/N C due to mixing of different operators hav- 
ing the same dimensions [44] . 

5. CONCLUSIONS 

• We have used the V-A component of the hadronic 
tau decays data for exploring the vacuum structure 
of the U LR QCD correlator using a set of Laplace 
sum rules (LSR). We have also revisited different es- 
timates based on FESR and its variant in section 4. 
Our results are summarized in Table 1. 

• Contrary to most papers in the literature, we do 



not perform a least square fitting procedure for ex- 
tracting simultaneously different condensates, but 
instead use the stability criteria (existence of minima 
or inflexion points) for our estimate of the conden- 
sates. Due to the alternate signs of the condensate 
contributions in the OPE and to the fact that in most 
methods, the high-dimension condensate contribu- 
tions are corrections to the lowest dimension conden- 
sates in the analysis, the approaches for extracting 
these high-dimension condensates can become inac- 
curate. 

• Instead, our strategy is to look for sum rules which 
disantangle, from the beginning, the relevant high- 
dimension condensates, and then makes the analysis 
cleaner and more transparent. 

• We have given a first estimate of the size of the 
D = 18 condensates, which will be interesting to 
check using alternative methods. The LSR estimate, 
which we expect to be more appropriate for extract- 
ing higher dimension condensates than wFESR and 
FESR, shows that the size of the very high D = 16 
and D — 18 condensates are relatively small which 
may indicate the good convergence of the OPE even 
at large r-values. 

• During the analysis, as one can see in previous 
figures, the absolute values of the condensates are 
slightly affected by r and t c in the optimum region 
(minimum or inflexion point). However, except for 
D = 4, it is important to notice that the results 
from LSR in large range of r- and £ c -values do not 
flip sign, which is a great advantage compared to the 
ones from some Finite Energy like-Sum Rules dis- 
cussed in the literature. 

• The extension of the present analysis to some 
other channels are feasible though not straightfor- 
ward. This is due to the relative importance of the 
continuum pQCD contribution for higher moments 
in some other channels, which is not the case of V-A, 
where this effect exactly cancels at higher energies. 
We plan to come back to these different channels in 
a future publication. 
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Table 1 

Estimated values of the D = 2d < 18-dimension (0_d) condensates in units of 1CP 3 GcV D . We have ordered 
the condensates from 0\& to O4 according to their chronological estimate. 



Authors 


018 


016 


014 


012 


010 


08 


06 


04 



THIS WORK 
LSR 

Eq. (8) -1± 0.6 +4.3 ±1.9 -9.6 ±3.1 +14.7 ±3.7 -17.1 ±4.4 +15.4 + 4.8 -9.7 + 4.1 -0.5 + 0.1 

Eqs. (14, 15) +15.8 + 3.2 -8.4+1.6 

Eq. (16) -8.0 + 1.1 

Average +15.6 ± 4.0 -8.7 ± 2.3 

wFESR 

Rev. 3 +30+ 10 -28 + 8 +25 + 5 -22 + 3 +16.8 + 2.0 -10.2 +.4 

Orig.[28] 4 -946 ± 147 +390 ±65 -U6±27 +43.5 ±10.5 -44 ±3.8 -4.8 ±0.9 

FESR 5 

BG Rev. 6 +12+1.5 -6.6 + 0.2 

BG Orig.[24] -12.4 ±9.0 -3.2 ±2.0 

DS Rev. 7 +10 + 2 -8.0 + 2.0 

DS Orig.[25] -2 ±12 -8.0 ±2.0 

LR Rev. 8 +53+ 16 -39 + 12 +26.0 + 8.4 -14.7 + 4.8 

LR Orig.[26] -260 ± 80 + 78 ± 24 -120+ 7 n -4 ±2 



OTHERS 



MHA+p'[19] +11.5 + 3.5-12.5 + 3.4+13.2 + 3.3-13.1 + 3.0 +11.7 + 2.6 -7.9+1.6 

MHA[19] +11.9 + 3.9-12.8 + 3.9+13.3 + 3.9-13.2 + 3.6+11.7 + 3.1 -7.9 + 2.0 

ZYA[22] -4.5 + 3.4 +7.8 + 3.0 -7.1 + 1.5 

IZ[23] +7.0 + 4.0 -6.4+1.6 

ALEPH[2] +11.0+1.0 -7.7 + 0.8 

OPAL[3] +7.5 + 1.3 -6.0 + 0.6 

DGHS[21] +8.7 ±2.4 -6.0 + 0.6 

CS3[29] -4.0 ±2.8 



AVERAGE 9 -1 ± 0.6 +14.4 ± 4.6 -15.7 ± 3.7 +23.8 ± 6.4 -18.2 ± 5.9 12.2 ± 2.9 -7.8 ± 1.6 



3 We have redone the analysis of [28] using instability criterion. 

4 We use the mean value of the results from the ALEPH and OPAL data. 

5 The revised (Rev.) FESR results have been obtained at t c 1.5 GeV 2 ; the original (Orig.) ones at t c fa 2.5 GeV 2 . 

6 These results have been obtained by [24] at 1.5 GeV 2 . 

7 We have corrected the value of Og (see section 4) and rescaled the results of [25]. 

8 The central values come from [27]. Inspired from the results of [26] at 2.5 GeV 2 , we have roughly estimated the 
systematic errors to be about 30%. 

9 Numbers in the lines Orig. are not considered into the average. 



